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Abstract 


Soil moisture is a crucial variable for weather prediction because of its 
influence on evaporation. It is of critical importance for drought and flood 
monitoring and prediction and for public health applications. The NASA Short- 
term Prediction Research and Transition Center (SPoRT) has implemented a 
new module in the NASA Land Information System (LIS) to assimilate 
observations from the ESA's Soil Moisture and Ocean Salinity (SMOS) satellite. 
SMOS Level 2 retrievals from the Microwave Imaging Radiometer using 
Aperture Synthesis (MIRAS) instrument are assimilated into the Noah LSM 
within LIS via an Ensemble Kalman Filter. The retrievals have a target volumetric 
accuracy of 4% at a resolution of 35-50 km. Parallel runs with and without 
SMOS assimilation are performed with precipitation forcing from intentionally 
degraded observations, and then validated against a model run using the best 
available precipitation data, as well as against selected station observations. 

The goal is to demonstrate how SMOS data assimilation can improve modeled 
soil states in the absence of dense rain gauge and radar networks. 


Description of LIS 



Description of SMOS & comparison to 
SMAP 



App l ications 

■ SPoRT produces real-time LIS soil moisture products for situational 
awareness and local numerical weather prediction over CONUS, 
Mesoamerica, and East Africa 

■ Currently interact/collaborate with operational partners on evaluation of 
soil moisture products 

o Drought/fire 
o Extreme heat 
o Convective initiation 
o Flood and water borne diseases 

■ Initial efforts to assimilate L2 soil moisture observations from SMOS 
(as a precursor for SMAP) have been successful 

■ Active/passive blended product from SMAP will be assimilated 
similarly and higher spatial resolution should improve on local-scale 
processes 


Description of Ensemble Kalman Filter 

We assimilate AMSR-E soil moisture observations using an Ensemble Kalman Filter (EnKF) 
within LIS. Kalman filtering is a data assimilation method that combines a forecast 
(background) with observations to generate an improved estimate of a model variable. A 
Kalman Filter calculates an optimal weighting between the background and the observation. 
The EnKF uses the spread of the ensemble to represent the forecast error covariance. We used 
an ensemble with 16 members generated using perturbations of 3 forcing variables (incident 
longwave and shortwave radiation, and rainfall), 4 state variables (4 layers of soil moisture), 
and 1 observation variable (SMOS soil moisture). 


Implementation of SMOS assimilation in LIS 

• Read ECMWF Level 2 Soil Moisture User Data Product (SMUDP) files 

• Read all orbits that fall in the time window (currently +/-3 hours) then exclude data 
outside the time window or the geographic region. 

• QC for RFI, frozen soil, snowcover, falling precipitation, heavy vegetation, and data 
quality flags 

• Use ESMF interface in LIS to assign observations to grid points. 

• Use LIS Ensemble Kalman Filter to produce model state updates. 





Experiment Setup 



"Truth" Run GDAS Single run Validation 
16 member ensembles. 

Spinup period of 8 months with state and forcing perturbations applied to 
generate ensemble spread. 

"Dry" and "Wet" runs with ±50% GDAS rain forcing are run in EnKF and 
open loop configurations to test the impact of SMOS soil moisture data 
assimilation on a known imperfect case as a proof of concept. 

Open loop runs are identical to EnKF runs with no data assimilated. 

Runs are validated against the run with the "correct" GDAS rain forcing. 
The EnKF runs assimilated data at 00Z and 12Z on 5 Jun 2013. 



Domain Eastern North America 

80x92 (25 km) 

Meteorological Forcing GDAS 

Model timestep 30 minutes 

Data assimilation timestep 6 hours (effectively 12 hours 

due to data availability) 


Background 

Model surface soil moisture 
(cm 3 /cm 3 ) at 12Z 5 Jun 2013 
before assimilation step. There 
was one prior data assimilation 
step at 00Z. 


Observations 

SMOS Level 2 retrieved surface * 
~12Z 5 Jun 2013. 


Case Study 


• Tropical Cyclone Andrea in June 2013 

• Before landfall, precipitation in Florida and Georgia associated 
with outlying rain bands 



GOES IR image. Courtesy of Naval 

Research Laboratory, Monterey, CA. 36-hr. antecedent 

precipitation from GDAS 
forcing data. 


Preliminary Model Results 



Discussion 

• This is a first test of SMOS assimilation so many shortcomings still need to 
be addressed. 

• Model "truth" (based on GDAS forcing and model physics) has errors which 
are not accounted for in this validation. 

• There is a large dry bias in the observations (relative to the model 
background) 

• Innovations are overall negative, even in the dry case, but there is a postive 
signal in the area of northern Florida/southern Georgia where precipitation 
occurred. 

• Analysis increments have the same sign as innovations, as expected. 

• Our relative error metric, hampered by model-observation bias, shows 
increasing error due to the analysis in most locations (red areas). However 
the total area with error reduction (green area) grows with repeated data 
assimilation cycles. 



Future Improvements 

• Apply a bias correction using the cumulative distribution function (CDF) adjustment of 
Reichle (???), possibly making it landcover-dependent (Blankenship and Crosson 
201 X) 

• Test assimilation in a cycling run over an extended period of time 

• Test in less-vegetated regions such as the Great Plains or tundra areas 

• Optimize the perturbations to achieve a spread that is representative of the estimated 

• Determine the optimal observation uncertainty 

• Perform validation against in situ soil moisture network observations. 

• Implement with SMAP test data available to SMAP Early Adopters 
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See also... 

■ Details on SPoRT's real-time LIS and expected use of SMAP data can be found in 

■ Zavodsky, B. T., J. L. Case, C. B. Blankenship, W. L. Crosson, K. D. White, 2013: 
Application of next-generation satellite data to a high-resolution, real-time land surface 
model, Earthzine, J. Kart, editor, Institute of Electrical and Electronics Engineers 
[Available online at http://www.earthzine.org/2013/Q4/10/application-of-next- 
generation-satellite-data-to-a-high-resolution-real-time-land-surface-model/ .l 

At this meeting: 

■ Case (26 th Conference on WAF/NWP; P162; Tues. Feb. 4) 

■ Case et al. (26 th Conference on WAF/NWP; J9.4; Wed. Feb 5 @ 2:15 P.M.) 




